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, In this paper, we study the chh'al symmetry restoration in the hadronic spectrum 

' in the framework of generahsed Nambu-Jona-Lasinio quark models with instanta- 

^ • neous confining quark kernels. We investigate a heavy-light quarkonium and derive 

o . 

, its bound-state equation in the form of a Schrodingerlike equation and, after the ex- 

cn 

' act inverse Foldy-Wouthuysen transformation, in the form of a Diraclike equation. 

o ■ 

' We discuss the Lorentz nature of confinement for such a system and demonstrate 

o 

fH , explicitly the effective chiral symmetry restoration for highly excited states in the 



mesonic spectrum. We give an estimate for the scale of this restoration. 
PACS numbers: 12.38.Aw, 12.39.Ki, 12.39.Pn 
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c5 ■ I. INTRODUCTION 

Chiral symmetry is spontaneously broken in QCD with light quarks. This immediately 
implies the existence of light mesons which become Goldstone bosons in the chiral limit, with 
dynamics described by chiral perturbation theory. Being an effective low-energy reduction 
of QCD, this theory does not describe the property of confinement. 

While the property of colour confinement is expected from general considerations and is 
supported by lattice calculations, no effective theory exists which derives confinement di- 
rectly from the QCD Lagrangian, and one is left to rely upon models, such as the constituent 
quark model, to describe confinement. These models work surprisingly well in the light- 
quark sector, reproducing the spectrum of low-lying mesonic states, with the exception of 
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ground-state pseudoscalars. So, the conventional wisdom of quark model practitioners is 
to assume that constituent quarks are indeed the correct degrees of freedom in the non- 
perturbative domain, with light pseudoscalars being a special case. As there is no hope to 
obtain Goldstone bosons within any naive quark model which does not address the question 
of QCD vacuum structure, the sector of light pseudoscalars has to be left outside the scope 
of the constituent picture. 

Phenomenological successes of constituent models are in fact quite remarkable, but there 
is yet another feature of QCD which they fail to reproduce. Namely, for highly excited 
hadronic states, we have quark typical momenta much larger than the scale of chiral sym- 
metry breaking so that chiral symmetry should be restored towards the high end of hadronic 
spectra. The phenomenon signalling this restoration is parity doubling. 

Parity doubling has been discussed for many years in connection with the baryon spec- 
trum (see, for example, jlf). Quite recently, the observation of the low-mass Ds(2317) and 
Ds(2460_]i mesons caused a resurgence of interest in the issue of chiral symmetry restora- 
tion 0, Detailed arguments in favour of parity doubling for higher mesonic states are 
given in pj. 

Quite obviously, chiral symmetry restoration is a consequence of fully relativistic dynam- 
ics of light quarks and chirally-symmetric quark interaction. Besides, it implies the confining 
quark interaction, so the models exhibiting parity doubling should interweave confinement 
and chiral symmetry breaking into one single mechanism. These models do exist in the 
literature flQ. and the possibUity of par.ty doubling in such models was briefly d,se„ssed 
in [21]. These models are introduced in Chapter II and generic mesonic bound-state equa- 
tions are derived for arbitrary kernels. Chapter III is devoted to the setting up of both the 
heavy-light bound state equation and its equivalent Diraclike equation. The heavy-light 
meson can be thought of as a kind of hydrogen atom of soft QCD, with the heavy-quark 
spin symmetry lOj, allowing us to study the dynamics of the light quark ignoring compli- 
cations related to the full relativistic two-body problem. Then, for arbitrary noncompact 
(that is, confining) quark kernels, it is shown in Chapter III that these equations necessarily 
lead to chiral restoration for high momenta. Finally, we end this chapter with a numerical 
study of the heavy-light mesonic spectrum, using, for the sake of clarity, a simple harmonic 
oscillator quark kernel (which allows for a differential mass-gap equation) and explicitly 
show the mechanism of chiral restoration at work. We conclude with Chapter IV. 
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II. MESONIC STATES IN POTENTIAL QUARK MODELS 



In this section, we briefly give the necessary details of the model to be considered below. 



The potential quar 



c model, which we call the Nambu.-Jpne-Lasinio (NJL)-type model after 



the original paper ll| is given by the Hamiltonian 



H 



d'x^ix, t)(-ii-S7 + m) ^(f , t) + ^ I d'xd'y J^x, t)Kf^{x - y)J'M t). 



with the quark currents, J"(x, t) = t)7^— t), interacting through an instantaneous 



kernel, 



Kl{x-y)=g,o9uoS'''Voi\x-y\) 



(2) 



A remarkable feature of the given model is the robustness of the results against variations 
of the form and parameters of the confining potential Vo{x). Usually a powerlike form is 
adopted. 



K.(|x|)--o 

with ^ a ^ 2. The case of the linearly rising potential a = 1 (see, for example 
strongly supported by hadronic phenomenology, whereas the marginal case of a = 2 



(3) 



13|) is 

— the 

harmonic oscillator potential — leads to simpler, differential, equations and is considered by 
many authors jo, LD, Isl llJ| because, despite of its mathematical simplicity, it already yields 
a physical picture for the mechanism of dynamical chiral symmetry breaking not unlike the 
one given by linear confinement. 

The quark models of this class fulfill the well-known low energy theorems of Gell-Mann, 
Oakes, and Renner [isl . Goldberger and Treiman Adler selfconsistency zero [l^]; the 
Weinberg theorem and so on. For an early derivation of the Gell-Mann-Oakes-Renner 
relation, for this class of models, see Ref. For the derivation of the Weinbergtheorem 
see Ref. Q]. The Salpeter equations for these models can be seen in detail in Ref. [7|. Using 
this formalism it is possible to give, for the class of models embodied in Eq. (P), an analytic 
proof of the Goldberger-Treiman relation 20]. The same formalism can be used to give an 
analytic proof of all the above low-energy theorems. The axial anomaly and the 7177 coupling 
constant can be derived naturally in the framework of the given model as well, following, 
for example, the lines of the textbook |2ll| and using the asymptotic freedom and the Ward 



identity for the dressed axial vertex 



(see also Ref. 



22| for an independent derivation of 



this relation in the Dyson-Schwinger approach to the dressed quark propagator). 
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The standard approach — vahd for any a — to the model (P) is the Bogohubov trans- 
formation from bare to dressed quarks, parametrised through the chiral angle y^p: 



-P)l (4) 




y/l + siuLfp + y/l - siny^p (aff) 
+ sinyjp - yT^^sirT^ (a^ 
where Ep stands for the dispersive law of the dressed quarks; (3 being the colour index, 



v{0), 



P = 1, Nc- In what follows the limit Nc ^ 1 is assumed. It is convenient to define the 
chiral angle varying in the range — | < ipp ^ ^, with the boundary conditions ip{0) = |, 
(p{p oo) 0. 

The Hamiltonian normal ordered in the basis (0)) can be split into three parts, 

H = E,^,+ :H2: + :H^:, (6) 

the first part being the vacuum energy, and the second and the third parts being quadratic 
and quartic in terms of single-quark operators. The requirement that : H2 '■ should be 
diagonal — the anomalous Bogohubov term be absent — leads one to the mass-gap equation 
for the chiral angle, 

Ap cos ipp — Bp sin = 0, (7) 

with 

Ap = m + ^ J -^;^V{p-k)sm(pk, Bp = p + ^ J {pk)V{p- k) cos ipk- (8) 

We absorb the fundamental Casimir operator Cp into the definition of the potential, V{p — 
k) = CpVo^p — k) by rescaling the parameter Kq. 

For the chiral angle — solution to the mass-gap equation ^ — the Hamiltonian ^ 
takes the form: 

^ = i^vac + E E / ^,Ep[%MhpM + dlX-p)dp,{-p)] + O {^-^=^ , (9) 

with the corrections coming from the : : part. This completes diagonalisation of the 
Hamiltonian (0) in the single-particle sector of the theory (the so-called BCS level 
The dressed-quark dispersive law can be built as 

Ep = Ap sin {pp + Bp cos ipp. (10) 



25[ in the axial 



It was found in Ref. l2Jl that the 't Hooft model for two-dimensional QCD 
gauge, which is identical to the model in two dimensions [26|, admits a further diagonali- 
sation in terms of colourless mesonic states, which is a step beyond the BCS approximation. 
This method was generalised, in Ref. [2^ , to the class of N JL-type four-dimensional models 
(^. The key idea was to rewrite the Hamiltonian (jHl), in the centre-of-mass frame. 



H = Kac + 

in terms of compound quark-antiquark operators, 

/ AT f ^ 



I3s\ 



(11) 



and to perform a second, generalised, Bogoliubov transformation, "dressing" the bare 
operators with the coherent cloud of quark-antiquark pairs. The operators creat- 
ing/annihilating the physical mesonic states in the theory are given by 





r p'^dp ' 


j rhnu = J 


(27r)3 [ 
f p^dp ' 




1 = j 


(27r)3 [ 



Mu{p)vtuip) - Ml{p)f-,ip) 

Ml{p)vt{p) ~ M,{p)VnAp) 

where the operators (|TT| are expanded first, 

MsAp^P) = Y}t^M]ss'Mu{p), 



(12) 



(13) 



using an appropriate basis diagonalising the spin-angular structure of the Hamiltonian. Such 
a basis is known to be formed by the {nJ^*^} states (n being the radial quantum number), 
for which we use the shorthand notation {nz/}. The centre-of-mass Hamiltonian li. takes a 
diagonal form. 



C 



(14) 



provided the mesonic wave functions obey the bound-state equation, 
[2E, - Mn,]^tAp) ' ^ 

with the orthonormality conditions 

p^dp 



i2.r 

k^dk 

T^r^siT-^P, k)vLik) + T--{p, k)^-^], 
(27r) 



(15) 



(27r)3 
p^dp 
(2^ 



[ptu{p)vLu{p) - ^nAP)^mv{P)\ = 
[ptnui.P)'^nui.P) - '^mui.P)VnAP)\ = 0. 



(16) 
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which follow directly from the canonical bosonic commutation relations of the operators ()12p. 
This completes the diagonalisation of the theory in terms of the physical observable 
degrees of freedom. The bound-state Ea. (ITHl) can be derived alternatively from the Bethe- 
S.pe.. e,„a.o„ fo. _ g ffl Q. We ... to .M. .et.od be.ow. T.e 

eigenvalue problem (fTHjl is subject to numerical investigation (see Ref. and also Ref. ^ 
where, in particular, the evolution of the vector-pseudoscalar mass splitting as a function 
of the current quark mass is studied in detail), and no problem is met when building the 
mesonic spectrum for the whole class of theories (^. 



III. A HEAVY LIGHT SYSTEM 



A. The bound state equation 



An important feature of the bound-state Eq. (fTH|) is the fact that each mesonic state 
is described by means of a two-component wave function, Lp^^{p). The physical meaning 
of these two components is obvious: ^Pnuip) describes the forward motion, in time, of the 
mesonic quark-antiquark pair and ^nuiv) time backward motion. Strictly speaking, for 
a given set of quantum numbers, four eigenfunctions should be considered — two for the 
eigenvalue Mny and two for —Mny. Meanwhile, the bound-state Eq. (fTH|l supports the 
symmetry 

{M„„ <^±,(p)} ^ {-M„„ vlM). (17) 

that can be easily verified using its explicit form. Due to the instantaneous interaction, 
both particles in the meson move forward and backward in time in unison. This leads us to 
an immediate conclusion that the heavy-light meson, with the static-quark Zitterbewegung 
suppressed by its infinite mass, can be described with a one-component wave function, so 
that 

Vnuiv) = V^nuip), '^nuip) = 0, (18) 

and the bound-state equation can be considerably simplified in this case. Although Eq. (fT3)) 
is written for the single-flavour theory, its generalisation for two flavours is trivial yielding 
for the heavy-light bound-state equation: 

/k'^dk 
-^T++{p,k)ipnuik), (19) 
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with niQ being the static-antiquark mass. The amphtude Tj"^{p, k) can be built using the 
Hamiltonian approach described in the previous section, though we find it more instructive 
to derive the heavy-hght bound-state Eq. (|19p using an alternative approach, employing 
the Bethe-Salpeter equation. To this end, we start from the homogeneous Bethe-Salpeter 
equation for the mesonic Salpeter amplitude x{P'i 



x{p;M) 



-Vip- k) ^oS.ik, ko + M/2)xik; M)SQ{k, k, - M/2)7o, 



where 



Sg{p,Po) 



A+(p)7o 



A-(p)7o 



Po- Ep + it po + Ep- ie 
is the dressed light-quark propagator. The projectors are defined as 



A±(p-)=TpP±Tt, Pi = 14^, T, = exp 



(20) 



(21) 



(22) 



Since, for the static particle, the chiral angle is simply ^q{p) = f , then the dispersive law 
and the projectors become trivial, giving for the heavy-quark propagator: 

P+lo , P-lo 



Sq{p,Po) 



+ 



Po-'mQ + te po + mq- te 
It is convenient to define a modified vertex. 



(23) 



r( 



^S,{p,po + M/2)x{p- M)Sq{p,po - M/2), 

ZTT 



(24) 



and to perform its Foldy-Wouthuysen transformation by means of the operator T^, 

r{p) = Tir{p)i, 



(25) 



where, for the sake of transparency, we kept the unity on the r.h.s. for the operator rotating 
the static antiquark. The resulting matrix bound-state equation reads: 

d^k 



(E - Ep)r{p) = 



(2vr) 



:V{p-k)TlTkr{k) 



(26) 



where we defined the energy excess over the static-particle mass, E = M — mq. The form 
of the r.h.s. of Eq. suggests that the matrix r{p) should have the structure 



r{p) 



i^ip) 




'ip{p) 



® (0 1)( 



(27) 
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where, for future convenience, we split the matrix into the direct product of the hght- and 
heavy-particle components. Then the eigenstate equation for the wave function tp{p) can 
be derived in the form of a Schrodingerlike equation. 



j^V{p- k) \C,C, + {ap){ak)S,Sk\ m = Eij{p), (28) 



with 



Cp = cos ^ - y^p) , Sp = sin ^ (I - <^p) ■ (29) 

The heavy-light bound-state Eq. is the main object for studies in the present paper. 
However, it is instructive to rewrite it in the Diraclike form, which allows one to investigate 
the Lorentz nature of confinement in the heavy-light quarkonium. First, notice that Eq. ()26|) 
is written for the positive-energy solutions, M = ttlq + E > 0. Similarly, for the negative- 
energy solutions, M = —mq — E < 0, although the same Schrodingerlike Eq. (j28|l holds, 
but the matrix r{p) takes a different form, 

m=i ° \ = f/M (30) 

which can be easily guessed from the symmetry p7j] and the specific property of the heavy- 
light system (fTHj). Thus the Foldy-rotated wave functions of the light component of the 
system, responsible for the positive- and negative-energy solutions, can be written, accord- 
ing to Eqs. (IHIUl . as 







for positive energies and as 





for negative energies. In view of the passive role played by the static particle, we shall 
refer to "^{p) as to the wave function of the entire heavy-light system. Consequently, the 
bound-state Eq. (j^Hj) can be rewritten for "^{p), 

Ep^{p) + J^^Vip- k) [TpTl + Tln\ m{k) = ^oE^{k). (31) 

As it follows immediately from Eq. (f^3j) . the Foldy operator which splits the positive- and 
negative-energy components of this wave function is T^: '^{p) = T^^{p). Therefore, the 
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pre-Foldy bound-state equation for the wave function '^{p) follows from Eq. ()31|) after the 
counter-rotation with the operator Tp, 

where, for future convenience, we introduced the unitary operator 

Up = Tp^oTl = 7o sin ipp + {dp) cos ^jp. (33) 

Then, using the mass-gap Eq. ((7j), rewritten in the form: 

EpUp = ap + l3m + ^J ^^(P - k)Uk, (34) 

we arrive at the heavy-light bound-state equation in the form of a Diraclike equation for 
the light quark, in coordinate space, 

{dp + /3m)^(f ) + j (fzk{x, z)K{x, z}^{z) = E^{x), (35) 

with 

K{x, z) = -[V^(f) + V{z) - V{x - z)] (36) 

and 

A{p,q) = i2n f5^'\p~q)Up. (37) 

^ ^ n 

Notice that the matrix A(p, k) |13i] can be related to the Green's function as 

A{p,k)=2t J ^S{u;,p,kho, (38) 
where we can use either the dressed light-quark propagator Sg, defined in Eq. (EH), 

S{u;,p,k) = {27Tf6^'\p- k)S,{p,u), 

or the entire heavy-light Green's function, built using the standard spectral decomposition, 
with the help of the solutions of the eigenstate problem (jH^ . 

E„>0 " E„<0 " 

This feature should not come as a surprise since, as mentioned before, the static antiquark 
decouples from the system, providing the overall colour neutrality of the bound state. 
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The Lorentz nature of confinement in the effective one-particle Eq. fl35l) follows from the 

, or equivalently, of Up for the given values of the relative 
momentum p involved. Below, we consider in detail the three regimes which can be identified 
in the spectrum of hadrons. Although the heavy-light system was used before to arrive at 
the quark-antiquark bound-state equation, all qualitative conclusions obviously hold for an 
arbitrary hadronic system. 

For small relative momenta the chiral regime takes place. Chiral symmetry breaking — 
spontaneous or explicit — plays a dominating role, the chiral angle being close to |. As a 
result. Up ~ 7o and the effective interquark interaction becomes purely scalar (even in the 
absence of any microscopic scalar force). For light quarks, the dressed-quark dispersive law 
Ep differs drastically from the free-particle form of ^/fP + m^, it may even become negative 
at small momenta — a necessary feature to have the light (massless) chiral pion. This 
chiral regime — where the effective interaction is purely scalar — holds below some effective 
dynamically generated low-energy scale, which we call the BCS scale Abcs, that is, holds 
for the mean relative interquark momentum p < Abcs- The BCS scale gives a measure of 
chiral symmetry breaking in the low-energy domain of the theory and it is closely related 
to the scale of confinement or, equivalently, to Aqcd^ Abcs — Aqcd- As soon as the mass 
of the quark-antiquark state grows, as does the mean relative interquark momentum, the 
vectorial part of the interquark interaction becomes more important. Finally, the chiral 
angle becomes negligibly small and the interaction acquires a purely vectorial nature. Thus, 
we arrive at the chiral symmetry restoration in the spectrum discussed in the literature 




Note, in passing, that the Dirac structure of the matrix A(p, k) is saturated by the scalar 
and space- vectorial part, so there is no room for the time-vectorial interaction, dangerous 
from the point of view of the Klein paradox. 



B. Chiral symmetry restoration in the spectrum 



1. General remarks 



It is argued in the literature [Jj, y, H ISl] that, for highly excited states in the hadronic 
spectrum, the chiral symmetry should be restored and the states with opposite parity should 
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come in doublets. This statement can be proved using the properties of the bound-state 
Eq. (j2Hl)- Indeed, let ip^p) be a solution to Eq. pHj) with the eigenvalue E. Consider 

^'(p) = m^ip), (39) 

which, by construction, possesses an opposite parity to iplp). The eigenvalue problem for 
the new wave function tp'{p) can be easily derived by multiplying Eq. ()28|1 by {ap) from the 
left, 

^\/(p - k) [S,S, + {ap){^k)C,Ck\ ij'ik) = Eij'ip), (40) 

where Cp and Sp are interchanged as compared to Eq. ()28|) . Notice that, for large momenta, 
one arrives at the relation Cp = Sp = which follows directly from the definition of the 
functions Cp and Sp fl29p and the large-momentum behaviour of the chiral angle, —>■ 
(see the first plot at Fig. 1). Therefore, for highly excited states, which possess a large 
relative interquark momentum, the bound-state Eqs. (j28j) and (I4(J|) coincide, being simply 

1 f d^k ^ r ^ ^ -1 ^ 

Epi^ip) + 2 y j^^^p- ^) + i^p)(^k)\ ^(^) = ^^(^^' (41) 

so that the states ipip) and {(ff))ip{p) become degenerate. Thus, in the framework of the 
potential quark model (PJ), one can investigate the explicit mechanism of the chiral symmetry 
restoration high in the hadronic spectrum. For further convenience, in addition to the BCS 
scale Abcs defined above, which can be related, for example, to the chiral condensate, 

3 f°° 

{qq) = 2 / dpp'^siuifp = -(ABcs)^ (42) 

^ Jo 

we introduce the scale of the symmetry restoration Aj-est? to be defined as the scale at which 
the splitting within the chiral doublet is much smaller than the BCS scale Abcs- Such a 
definition is natural since, alternatively, the BCS scale can be defined as the splitting within 
a chiral doublet at small values of the relative interquark momentum, that is, in the region 
where this splitting is essentially dictated by chiral symmetry. The numerical estimate for 
Arcst in the case of the harmonic confinement {a = 2) is given below. 

For large interquark momenta, the functions Cp and Sp are smooth and slow, remaining 
almost constant when approaching their large-p asymptote of Therefore, since, for 
confining potentials, the distribution given by the Fourier transform V{p—k) has the support 
at /c ~ p, it is possible to consider, approximately, CpCk ~ Cp, SpSk ~ 5"^ and to study Cp 
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1 0,0 0,5 1,0 1,5 2,0 



P P 

FIG. 1: Solutions to the mass-gap Eq. (O (the first plot) and the coefficients and Sp (the 
second plot) for the potential V{r) = K^'^^r"' with various a's. We plot the curves for a = 0.3, 
0.5, 0.7, 0.9, 1.0, 1.1, 1.3, 1.7, and 2.0 (smaller a's correspond to steeper curves at the origin). The 
momentum p is given in units of Kq. 

and Sp (see the second plot at Fig. 1). This approximation can be easily checked for the 
harmonic oscillator potential, for which V{p— k) oc A^d^^^p— k). Taking the corresponding 
integral by parts, one can omit, for large momenta, the terms containing derivatives of 
the almost constant functions Ck and Sk- Thus, the chiral symmetry is restored as soon 
as the difference Cp — Sp becomes negligible. Using the definition one can find that 
Cp — Sp = simpp = f^{p), with f^{p) being the wave function of the chiral pion, with 
coinciding positive- and negative-energy components 0, Q, Isl- Therefore, we conclude 
that chiral symmetry is restored in the spectrum if the pionic wave function vanishes for 
the given values of the relative interquark momentum. This gives us an explicit physically 
transparent relation between the BCS scale and the scale of the symmetry restoration. 
Although, formally, both scales Arest and Abcs are defined by the only dimensional parameter 
in the theory Kq, we expect, due to numerical factors, the relation Ajcst ^ Abcs to hold. 
Below, we perform a numerical check of this relation and estimate the value of the restoration 
scale. 
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2. Numerical estimates: Harmonic oscillator potential 

In this paragraph and in order to illustrate the general results obtained so far, we per- 
form an explicit numerical study of the heavy-light mesonic spectrum. To avoid unnecessary 
complications, we consider the case of the harmonic oscillator potential, which allows us to 
formulate the mass-gap equation and the bound-state problem through differential equa- 
tions. This corresponds to the case a = 2 for the potential Q. Then the Fourier transform 
of the potential is the Laplacian of the three-dimensional delta-function, 

V{p- k) = -K^Ak6'~'\p~ k), (43) 

the chiral angle is the solution to the differential mass-gap equation, 

sin (Pp = ^Kq [p'^ipp + 2pip'p + sin 2(pp\ + mp^ cos v^p, (44) 

and the dressed quark dispersive law can be calculated as 



. 3 '{V'p? , COsVp 

tjp = m sm ^p + p cos (pp — ' 



+ 



(45) 



2 ■ p2 

In order to introduce the BCS scale we, following the approach suggested before, evaluate 
the chiral condensate, 

{m) = i / dpp'sin^p ^ -(O.SliTo)' = -(Abcs)', (46) 

^ Jo 

which gives us 

= mMeV, (47) 

if the standard value of the chiral condensate of —{250MeV)^ is used in order to fix Abcs- 
Following the definition of the symmetry restoration formulated above, we demand that the 
mass splitting in a chiral doublet AE should be much less than Abcs? being, at most, just 
a few per cent of 250MeV. 

In order to solve the bound-state Eq. ()28|1 . we use, as a conventional scheme for the 
quantum numbers u, the angular momentum I and the total momentum j = / ± | of the 
light quark, so that the wave function ip{p) can be expanded as 

m = ^jiM^^ (48) 
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with u{p) being the radial wave function, and the spherical spinors are defined as 



^Jlmif^ = 5Z^/wiM2^'^l(^^^- 



(49) 



Then, the eigenvalue equation for the radial wave function u{p) can be readily derived 
from the bound-state Eq. (j28p in the form: 



u 



" = [Ep - E]u + Kl 



K 

+ + ^ sm ipp 



(50) 



where 



/, for j = / - i 
-(/ + 1), for j = / + i 



Notice that Eq. (|Kn|l can be rewritten in the form of a Schrodigerlike equation. 



(51) 



with the effective potential 



,2 , (j + l/2f K . 

— + — sm ifp 

pZ pZ 



(52) 



The well-known property of the spherical spinors. 



{a^p)Vtjim{Pj = -^Ijl'mif), l + l' = 2j, 



(53) 



ensures that the states with j = I ± ^ possess opposite parity. The splitting between such 
states is the subject of our investigation, and it follows immediately from Eq. (jK^ that, for 
the given total momentum j, this splitting is due to the ^-dependent term in this effective 
potential. Thus, for k = ib(j + |), one finds for the difference of the opposite-parity 
potentials: 

-J- 1 

(54) 



p. 

which gives us an explicit example of the relation between the pionic wave function sin^^p 
and the opposite-parity state splitting, discussed above in general terms. 

It is clear from Eq. (j54p and from the form of the chiral angle (the first plot at Fig. 1) 
that the splitting vanishes for excited states, since the excited wave function is localised at 
larger relative momenta, whereas the chiral angle decreases fast with p. In order to prove 
this property explicitly, we solve the bound-state Eq. numerically. The results are listed 
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TABLE I: The masses of orbitally excited states and the sphttings for the radial quantum number 
n = 0, as solutions to the bound-state Eq. H51() with the potential (|52() . The solutions of the 
Salpeter Eq. H55|) are also listed for the sake of comparison. All energies are given in units of Kq. 

j 1/2 3/2 5/2 7/2 

E,_.i 2.04 3.51 4.51 5.35 

E,_,.i 2.66 3.69 4.57 5.36 

AEj 0.62 0.18 0.06 0.01 

1 2.34 3.36 4.24 5.05 

^=3-2 

Ef'^'P 1 3.36 4.24 5.05 5.79 

AE^""^"^ 1.02 0.88 0.81 0.74 



TABLE II: The masses of orbitally excited states and the splittings for the radial quantum number 
n = 1, as solutions to the bound-state Eq. with the potential The solutions of the 

Salpeter Eq. H55|) are also listed for the sake of comparison. All energies are given in units of Kq. 

j 1/2 3/2 5/2 7/2 

E,_-i 3.91 5.03 5.87 6.60 

E,_.,i 4.39 5.17 5.92 6.61 

'— J^2 

AEj 0.48 0.14 0.05 0.01 

-Ef^'P 1 4.09 4.88 5.63 6.33 

Ef'""^^ 1 4.88 5.63 6.33 7.00 

A£;|^'P 0.79 0.75 0.70 0.67 



in Tables I,n. They demonstrate a clear pattern of chiral symmetry restoration for orbitally 
excited heavy-light mesons. For the sake of comparison, we calculate also the spectrum of 
the Salpeter equation 

[^/p^ + m? + Klx^]i){x) = Ei){x), (55) 

which, in momentum space, also has the form of Eq. (j5H) with the potential (|52p. with the 
quark dispersive law substituted by the free energy and the chiral angle put equal to | in 
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FIG. 2: Regge trajectories for Eq. H51() with the potential (|52p (soHd hnes) and for the Salpeter 
Eq. (dashed hnes). In each pair, the lower curve corresponds to / = j — ^ and the upper curve 
corresponds to I = j + ^. 

the interaction, which is, 

v|f;;p(p) = V^T^ + Ki ^^^^'^f^" = V^T^ + K^-^^. (56) 

The results are also listed in Tables I, II. 

In Fig. 2 we plot the Regge trajectories for the bound-state Eqs. and as E^^'^ 
versus j^, for the radial quantum numbers n = (the first plot) and n = 1 (the second 
plot). For the bound-state Eq. (jKT|) . the trajectories for j = / ± | merge in such a way that 
at j = 5/2 the splitting is 0.06 (in the units of Kq), whereas for j = 7/2 it is already 0.01. 
Now the BCS scale Abcs is 0.51Kq. Therefore, if we require that the splitting constitutes 
a few per cent of 0.51 (say, about 10 per cent), then the splitting should be below 0.051, 
which happens between j = 5/2 and j = 7/2, at E^^^ ^ 11. This gives the bound-state 
energy E ^ 5Kq. Let us define this energy to be the restoration scale Arest- For j = 7/2 
we have a mass splitting which, despite being already fifty times smaller than Abcs, yields 
a bound-state energy (of ~ 5.3) just slightly above Arest- It is then clear that Arest can be 
thought of as the onset scale for the chiral restoration regime. In other words: for a mass 

^The quasiclassical spectrum of Eq. (|55|l gives E oc p/^. Therefore, for this equation, for the given radial 
excitation number n and the total momentum j, one expects two neighbouring nearly parallel straight line 
trajectories E^^^{j), corresponding to / = j ± i. Fig. 2 clearly exhibits such a behaviour. 
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FIG. 3: The ratio of the sphtting Ai?,- = E,_-,i — E,_-i to the BCS scale AbcSj in per cent, 
versus the averaged bound-state mass Ej = ^ (^Ei^j^i + Ej^^j_i^. For the sake of transparency 
we measure the bound-state energy in GeV, using the estimate (|17)) . 

splitting 10 times smaller than the BCS scale Abcs we have A^cst 10 times bigger than the 
BCS scale Abcs; 

Arest ~ 5Ko ^ 2.5GeV ^ IOAbcs, (57) 

where the estimate ipTj) was used for the parameter Kq. The result fl57j) is in good agreement 
with other estimates of the restoration scale known in the literature In Fig- 3 we 
illustrate, for the sake of completeness, the behaviour of the splitting AEj = Ei^-j^^_—Ei^-_i 
as a function of the averaged bound-state energy Ej = ^{Ei^j_^_i + Ei^j_i). 

In contrast, the trajectories for Eq. remain parallel for all values of the total momen- 
tum j. This contrast provides an explicit demonstration of the role played for the restoration 
of the symmetry by the pionic wave function siiupp in the potential (jK^ . 

3. Generalisation to a generic form of the potential 

The results obtained in the framework of the potential quark model are known to be 
robust against variations of the form of the quark kernel — although the values of various 
physical quantities might change, for example, with the change of the shape of the confining 
potential, the relations between such quantities are left intact. 
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FIG. 4: The coefficients and Sl for the potential V{r) = K^^+^r'^ with various a's (we use the 
same set of a's as in Fig. 1). 

In Fig. 4, for the sake of clarity, we plot the coefficients Cp and 5*^ versus the relative 
momentum p, measured in GeV. To this end we extracted the corresponding scale Kq, for 
the given power a, from the definition of the BCS scale — see Eqs. (j46|) . (|47|). This plot is 
to be compared with the second plot in Fig. 1, where the same quantities are plotted in the 
units of Kq. One can see the universality of the behaviour of these functions, regardless of 
the explicit from of the confining potential. A similar conclusion holds for the dressed quark 
dispersive law Ep which, for excited states in the spectrum, tends to the a-independent free- 
quark limit of a/p^ + m^. We conclude, therefore, that the only ingredient of the bound-state 
Eq. ()28|1 which changes with the change of the confining potential is the Fourier transform 
V{p — k), the net result of this change being simply the modification of the Regge trajectory 
form, from E^^'^ vs j, for the harmonic oscillator potential, for E^ vs j, for the generic 
potential 0. Therefore, we dare to extend the relation (jFTjl beyond the harmonic oscillator 
potential case of a = 2, generalising it to the entire family of allowed confining potentials. We 
find it remarkable that the chiral potential model suggested for studies of chiral low-energy 
phenomena in QCD appears able to address the problem of the chiral symmetry restoration 
in the spectrum of highly excited hadrons and to give a reasonable robust prediction for the 
corresponding restoration scale. 
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IV. CONCLUSIONS 



In this paper, we investigate the chiral symmetry restoration for highly excited states in 
the hadronic spectrum. We consider the NJL-type model with an arbitrary confining 
quark kernel, and use it to study the heavy-light quarkonium. We derive the bound-state 
equation for the given system in the form of a Schrodingerlike equation and also in the form of 
an effective Diraclike equation for the light quark in the field of the static antiquark. We give 
an explicit expression for the matrix structure of the effective interaction of the light quark 
with the antiquark source for all relative interquark momenta and study the Lorentz nature of 
this interaction. Thus, we identify explicitly the following three regimes: i) the chiral regime, 
P ^ ^BCS) with chiral symmetry breaking playing a dominating role in the interaction — 
the latter being predominantly scalar, ii) the restoration regime, Aj-est ^ P, which realises 
for highly excited states in the spectrum, with the effective interaction being vectorial and, 
as a result, states with opposite parity coming in doublets, and iii) the intermediate regime, 
which interpolates the two aforementioned ones. We perform a numerical investigation of 
the heavy-light bound-state equation for the harmonic oscillator potential and demonstrate 
explicitly the effect of chiral symmetry restoration in the spectrum for orbitally excited 
states. We estimate the restoration scale to be Arest ~ IOAbcs ~ 2.5GeV, where the BCS 
scale AbcS) defining the low-energy properties of the theory, for example, related to the 
chiral condensate, is chosen to take the standard value of about 250MeV. Bearing in mind 
robustness of predictions of the theory for relations between various physical parameters 
with respect to variations of the quark kernel, we dare to extend this conclusion to an 
arbitrary confining interquark kernel in the Hamiltonian (P). Thus, we conclude that there 
is a sufficient window for the intermediate regime in which, on one hand, chiral symmetry 
does not play a dominating role anymore, whereas, on the other hand, the parity doubling 
still does not happen in the spectrum. We find that the scale of the chiral symmetry 
restoration in the spectrum, evaluated in the framework of the potential models (^, is in 
good agreement with other estimates known in the literature. 
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